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Abstract. This paper addresses the theoretical analysis of dilepton spectra in 
C+C collisions at GSI-SIS energies. Theoretical predictions for the recent data 
of the HADES collaboration at SIS-energies are made with the help of a hadron- 
string transport model, the Ultra-relativistic Quantum Molecular Dynamics 
(UrQMD) model. A mass shift of the p meson due to kinematical effects is 
discussed. 
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1. Introduction 

Dilepton spectra are expected to play an important role in indicating deconfincment 
and the restoration of the spontaneous breaking of chiral symmetry. Particularly 
in the low mass region, dileptons couple directly to light vector mesons and can 
therefore reflect the mesons mass distribution in the hot and dense region. Since 
dileptons do not interact strongly with the surrounding hadronic medium, they 
escape the fireball of a heavy ion collision and might carry information on properties 
of the early stage of the collision. Thus, they are one of the most crucial observables 
for the exploration of the QCD mass generation [ ^ |2J |3J 0J |S] • 

Heavy ion reactions with energies around 1 — 2 AGeV create a reaction zone 
of high nuclear density for about 10 fm/c. This is long compared to the typical 
lifetime of a hadronic resonance. Thus, the p-mesons with lifetimes of r ~ 1 fm/c, 
can decay within this zone and the decay products (i.e., the e + e~-pairs) still carry 
the properties of the dense medium. The u> meson (r ~ 22 fm/c) and the <j) meson 
(r ~ 44 fm/c) have longer lifetimes and can therefore escape the fireball before 
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they decay. Consequently, the p meson is the most prominent candidate to study 
in-medium modifications of hadron properties in high density nuclear collisions. 

Calculations at low beam energies fall in the category of non-perturbative QCD, 
i.e., no exact methods to solve the equations of motion at QCD exist. Therefore 
one is forced to utilise approximations within effective field theories and compare 
with the results of experiments. In this paper, we focus on the HADES experiment 
at GSI-SIS [E|. 

For our studies we apply the UrQMD model. It is a non-equilibrium transport 
approach based on the covariant propagation of hadrons and strings. All cross 
sections are calculated by the principle of detailed balance or are fitted to available 
data. The model allows to study the full space time evolution for all hadrons, 
resonances and their decay products. This permits to explore the emission patterns 
of the resonances in detail and to gain insight into the origin of the resonances. 
UrQMD has been successfully applied to study light and heavy ion reactions at 
SIS. Detailed comparisons of UrQMD with a large body of experimental data at 
SIS energies can be found in [Cj. For further details of the model the reader is 
referred to [EB Ej. 



2. Theoretical Background for Dilepton Production 

2.1. Dalitz Decay of n°, rj, u and r{ 

Here, we follow [E3] and [ lllj to derive the equations for the decay rate of a meson 
into dileptons. In general, these processes have the form 

P^Ve+e-, V^Pe+e-, (1) 

with P being a pseudoscalar meson and V being a vector meson. There processes are 
not true three-body decays because the processes can be split into two subprocesses. 
Thus, the decay A — > Be + e~ can be separated into the decay of a virtual photon, 
A — > B7*, which subsequently decays via electromagnetic conversion, 7* — ► e + e~. 
According to [|H3| and [E] we arrive at the Dalitz decay formulas, which differ only 
in the values for the form factors and if the particle B is a real photon (to b — 0) or 
a virtual photon (mg 7^ 0): 

• Dalitz decay of the pseudoscalar mesons ir°, rj and rf (tob = 0): 



dM 3irM V M 2 \ M 2 



dN A _^ e+e - _ 4a 4m 2 ( 2m 2 \ { ^ M 
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x \F AB (M 2 )\ 2 ^^(N A ). (2) 

1 tot 

The quantity (Na) is the number of mesons A per event, M is the invariant 
mass of the lepton pair, m e the electron mass, and F AB the corresponding 
form factor. 
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Dalitz decay of the vector meson oj (tub 7^ 0): 
dM 




(Na) 
(3) 



The Dalitz decay of the A resonance, A — » Ne + e~, differs from the previous 
discussed decays because of the different interaction NAj vertex. The decay rate 
reads as follows 

dN e+e - _ f dN A ^ Ne+e - rn f s dN Ajlljr _ f 2a T(M A ,M) dN A 

->tot 
AO 



(M A ) — —dM A — / y -^- f — '- — —dM A , (4) 

dM J dM y AJ dM A J 3ttM T f fn dM A ' w 



where the fraction dN& 2u + "~ ( M a) = ^tT^f ( M <\) is the number of e+e~- 
pairs per A resonance, invariant mass bin and event. M A is the actual mass of the 
A resonance from the simulation and T A q is the total decay width at the resonance 
pole mass. Additionally, we use the decay width into a massive (virtual) photon 

T(M A ,M) = L ' f A) m N [2M t (M, M A ) + Mi(M, M A )]. (5) 

1ottM a 

The function A is defined by A (m A , mf, m|) = (m A — (mi+m 2 ) 2 )(m^- (mi — m2) 2 ). 
The corresponding matrix elements .Mt and A^; are taken from [I12j. 

2.2. Direct Decay 

The decay of a neutral vector meson in an e + e~-pair is a true two-body decay. 
Thus, energy and momentum conservation leads to the fact that the mass of the 
original meson (and the virtual photon) is the same as the invariant mass of the 
dilepton (my = M e + e ). For a direct decay the branching ratio reads [ 11 BR(V — > 
e + e~) = r v^+e- ( M ) with V = p°, u>, (f>. T v _, e + e - varies with the dilepton mass 
like M~ 3 according to [ I13| 



T v ^ e+e - (A/) — — y 1 - — ^1 + 2— j (6) 

with T v ^ e + e - (my) the partial decay width at the resonance peak mass. For each 
meson mass the branching ratio (BR) is calculated and summed over in each mass 
bin AM to get the distributions dN/dM for the figures shown below. 

The present calculations do not consider contributions from bremsstrahlung, 
because of the relatively small rates [114]. 
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2.3. Implementation into the Simulation 

In order to extract the dilepton abundances from the UrQMD model [ [21 El, we 
use the lifetimes and 4-momenta of the calculation for each channel. A meson can 
be produced via the decay of a meson or baryon resonance, annihilation or string 
fragmentation. Such a meson obtains its mass according to the relativistic Brcit- 
Wigner distribution, and a 4-momentum vector resulting from the kinematics in 
the single process under consideration. After the creation, the meson is assumed to 
travel on a straight line until it decays or it collides with another particle and its 
4-momentum changes. 

The most time consuming factor with the numerical analysis of electromagnetic 
processes is the suppression by a factor of a 2 « 10 -5 due to the electromagnetic 
coupling. Therefore, we assume that the vector mesons radiate dileptons continu- 
ously during their whole lifetime ("shining" method). After the global freeze-out 
time tf, the dilepton yield can be treated as in [El and [E] an d we get 



The second term of Eq. JJJ) describes the process very well if the decay rates do not 
depend on time. The invariant mass distribution of the dileptons is then equal to 
the distribution of the meson (V) itself multiplied with the branching ratio. But in 
a time dependent simulation one has to consider collisional broadening, i.e., the fact 
that the dilepton can be reabsorbed before it decays. This means that the lifetime 
t of the vector meson producing the lepton pair and the decay width T varies for 
each considered hadron. In this case it is more suitable to use the time integration 
method. 

3. Dilepton Yields at SIS Energies 

In this section we present predictions for dilepton spectra which correspond to the 
experimental runs carried out by the HADES collaboration [El- We consider mini- 
mum bias (b < 4.5 fm) C+C collisions at 2 AGeV. For the calculation we utilise the 
UrQMD model. To achieve better statistics, we randomly create 1 to 100 dilepton 
pairs per hadron decay (depending on hadron species and original abundances) with 
the direction and momenta of the dileptons in the local rest frame of each hadron. 
These e + e~ momenta are transformed to the laboratory frame. For convenience, 
all spectra are then normalised to the mean pion multiplicity (N^o). 

The left hand side of Fig. ^depicts the invariant dilepton mass spectra of C+C 
collisions at 2 AGeV. The pion multiplicity obtained from the UrQMD calculation 
is (N^o) — 1.223/event. One observes that for masses below 110 MeV the 7r°-Dalitz 
decay contribution is up to two magnitudes higher than the other contributions but 
does not account for higher masses. The A- and 77-Dalitz decay outweigh the other 
yields until 400 MeV. For masses from 400 to 600 MeV the w-Dalitz and direct p 



dN e+e 
dM 




(M) + 



Y e+e -{M)dN v {t f ) 



(7) 
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decay contributions become important. The direct decays of the p and also that of 
the u) are responsible for the yields with the highest masses up to 1 GeV. 
We include the coupling of the p meson to pions and baryons via the employed 





Fig. 1. Invariant mass spectra for dilepton abundances in minimum bias C+C 
collisions at 2 AGeV (left). The yields are normalised to the pion multiplicity (N^o) 
= 1.223/event. In the right panel the p meson is distributed via a free Breit-Wigner 
distribution. 

cross sections fitted to available data or calculated from detailed balance [0 E|- In 
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the intermediate-mass region around 500 to 600 MeV the p spectrum is enhanced 
due to the decay chain 7V*(1520) — > N + p; p — > e + e~. For higher p masses heavier 
baryon resonance decays and pion-pion annihilation are the most contributing parts 
to the spectrum. This leads to a small second hump in the dilepton spectrum. One 
hump is centered around 500 to 600 MeV and one at the p meson vacuum mass 
(m p = 770 MeV). 

Although the present calculation is based on vacuum cross sections and vacuum 
hadron widths the spectra from this many body dynamics look like what is expected 
from calculations employing in-medium spectral functions explicitly Due to the 
decay kinematics of the production channel of the p meson a strong modification 
of the p spectral function occurs. This effect is of kinematic origin only and will be 
discussed in the following section. 

To investigate the difference between the p meson distribution originating from 
the UrQMD simulation with all the kinematics applied and the vacuum one, we 
compare now the spectra above to the case that the p mass has a free Breit-Wigner 
distribution. Figure^ (right) depicts this scenario, the other rates remain unaltered. 
The discrepancy is obvious: The low-mass part of the p is suppressed whereas masses 
over 900 MeV appear. Additionally, a clear maximum around the p vacuum mass of 
770 MeV is visible. This leads to a dip in the sum of all rates around 600 MeV. The 
original UrQMD spectra, in contrast, show no maximum at the vacuum mass. It is 
an approximately flat curve for 2 AGeV. In conclusion, the present transport model 
simulations suggest strong kinematic effects and leads in fact to modifications of 
the p mass spectrum and dilepton spectra. 

Fig. |2|depicts invariant dilepton mass spectra from C+C collisions at 2 AGeV cor- 
rected to the HADES acceptance. The calculated dilepton momentum vectors were 
filtered by the HADES collaboration (the official version of the HADES acceptance 
filter, release 1.0, was used). This filter contains an acceptance matrix for each 
point at full solid angle and a smearing of the momenta according to experimental 
resolution. 



4. p meson contributions 

Let us finally discuss the p meson contributions to the dilepton spectrum. Fig. |3do- 
picts the p meson mass spectrum for minimum bias C+C reactions at 2 AGeV beam 
energy. One observes a double peak structure, where the peak at ~ 500 MeV origi- 
nates from A*(1520) resonance decays, i.e., by the decay chain AT*(1520) — > N + p. 
As seen in Fig. p mesons originating from such a baryon resonance decay have 
masses below 600 MeV. This effect is of kinematic nature only, since the calculation 
does not make use of explicit in-medium spectral functions, p mesons originating 
from 7T7T annihilation or higher mass baryon resonance decays have a mass, which is 
compatible with the pole mass of 770 MeV. Note that the p production is dominated 
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Fig. 2. Invariant mass spectrum for dilepton production in minimum bias C+C 
collisions at 2 AGeV. The yields are filtered by the acceptance matrix of the HADES 
collaboration, release f .0. 

by the decay of baryonic resonances. Direct p production from ttw scattering yields 
only a small contribution up to 20 % at the p-peak mass. For more details regarding 
p mesons in C+C collisions at 2 AGeV beam energy, the reader is referred to [I17|. 

5. Conclusions 

We have shown that due to the collision kinematics the mass spectrum of vector 
mesons, especially of the p, is modified. Therefore, a many body transport approach 
without explicit use of in-medium spectral functions leads to sizeable shifts in the 
dilepton mass distribution compared to a vacuum baseline. The prediction of a 
kinematically induced mass shift at 2 ^4GeV beam energy can be tested by the 
HADES experiment, e.g. when going towards target/projectile rapidities. 
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